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Pure hadronic compact stars, above a threshold value of their gravitational mass (central pressure), 
are metastable to the conversion to quark stars (hybrid or strange stars). In this paper, we present 
a systematic study of the metastability of pure hadronic compact stars using different relativistic 
models for the equation of state (EoS). In particular, we compare results for the quark-meson 
coupling (QMC) model with those for the Glendenning-Moszkowski parametrization of the non- 
linear Walecka model (NLWM). For QMC model, we find large values (M cr = 1.6 - 1.9M Q ) for the 
critical mass of the hadronic star sequence and we find that the formation of a quark star is only 
possible with a soft quark matter EoS. For the Glendenning-Moszkowski parametrization of the 
NLWM, we explore the effect of different hyperon couplings on the critical mass and on the stellar 
conversion energy. We find that increasing the value of the hyperon coupling constants shifts the 
bulk transition point for quark deconfinement to higher densities, increases the stellar metastability 
threshold mass and the value of the critical mass, and thus makes the formation of quark stars 
less likely. For the largest values of the hyperon couplings we find a critical mass which may be as 
high as 1.9 - 2.1 Mq. These stellar configurations, which contain a large central hyperon fraction 
{fy,cr ~ 30%), would be able to describe highly- massive compact stars, such as the one associated 
to the millisecond pulsars PSR B1516+02B with a mass M = U^^Mq. 



PACS number(s): 97.60.S, 97.60.Jd, 26.60.Dd, 26.60.Kp 



I. INTRODUCTION 

The nucleation of quark matter in neutron stars has been studied by many authors, due to its potential connection 
with explosive astrophysical events such as supernovae and gamma ray burst. Some of the earlier studies on quark 
matter nucleation (see e.g., [l], 0, H| and references therein) dealt with thermal nucleation in hot and dense hadronic 
matter. In these studies, it was found that the prompt formation of a critical size drop of quark matter via thermal 
activation is possible above a temperature of about 2 — 3 MeV. As a consequence, it was inferred that pure hadronic 
stars are converted to quark stars (hybrid stars (HyS) or strange stars (SS)) within the first seconds after their 
birth. However, neutrino trapping in the protoneutron star phase strongly precludes the formation of a quark phase 
0) IS B 13 • Then, it is possible that the compact star survives the early stages of its evolution as a pure hadronic 
star. In this case, the nucleation of quark matter would be triggered by quantum fluctuations in degenerate (T = 0) 
neutrino free hadronic matter 

H ©, 53, ED, E El, 03, El- 
Quantum fluctuations could form, in principle, a drop of /3-stable quark matter (hereafter the Q 13 phase). However, 

this process is strongly suppressed with respect to the formation of a non /3-stable drop by a factor ~ Cj^rmi , where 
N ~ 100 — 1000 is the number of quarks in a critical-size quark drop. This is so because the formation of a /3-stable 
drop would involve the almost simultaneous conversion of ~ A/3 up (u) and down (d) quarks into strange (s) quarks. 
Alternatively, quantum fluctuations can form a non /3-stable drop (hereafter the Q* phase), in which the flavor content 
of the quark phase is equal to that of the /3-stable hadronic phase at the same pressure [3, [13, G3 • Since no flavor 
conversion is involved, there are no suppressing Fermi factors, and a Q* drop can be nucleated much more easily. 
Once a critical-size Q* drop is formed, the weak interactions will have enough time to act, changing the quark flavor 
fraction of the deconfined droplet to lower its energy, and a drop of the Q" phase is formed. This first seed of quark 
matter will trigger the conversion 

[13, 53, EDI of the pure hadronic star to a hybrid star or to a strange star (depending 
on the details of the equation of state for quark matter used to model the phase transition). The stellar conversion 
process liberates a total energy of the order of 10 53 erg [l8| . 

When finite-size effects at the interface between the quark and hadron phases are taken into account, it is necessary 
to have an overpressure AP = P — P Q > with respect to the bulk transition point Pq, to create a drop of deconfined 
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quark matter. As a consequence, pure hadronic stars with values of the central pressure larger than Pq are metastablc 
to the decay (conversion) to hybrid stars or to strange stars [Til . [l2l . [T3L EH, EH] • The mean lifetime of the metastable 
stellar configuration is related to the time needed to nucleate the first drop of quark matter in the stellar center and 
depends dramatically on the value of the stellar central pressure [n], EU, [13, EU, EH • 

The possibility of having in nature both metastable hadronic stars and stable quark stars, has led the authors of 
ref. [l2| to extend the co ncep t of limiting mass of a "neutron star" with respect to the classical one introduced by 
Oppenheimer and Volkoff [19| , Since metastable HS with a "short" mean-life time are very unlikely to be observed, 
the extended concept of limiting mass has been introduced in view of the comparison with the values of the mass of 
compact stars deduced from direct astrophysical observation (see sect. 3.1 of ref. 12] for the definition of the limiting 
mass, Mu m , of compact stars in the case of metastable pure hadronic stars). 

As it is well known, neutron star mass measurements give one of the most stringent test on the overall stiffness of 
dense matter EoS. Recent measurements of Post Keplerian orbital parameters in relativistic binary stellar systems 
(containing millisecond pulsars) give strong evidence for the existence of highly-massive "neutron stars". For example, 
the compact star associated to the millisecond pulsar PSR B1516+02B in the Globular Cluster NGC 5904 (M5) has 
a mass M = 1.94±q \IMq (1 a) [20] . In the case of PSR J1748-2021B, a millisecond pulsar in the Globular Cluster 
NGC 6440, the measured mass is M = 2.74~^q'^\Mq (2 a) [2l[. These measurements challenge most of the existing 
models for dense matter EoS. 

In this work, we carry out a systematic study of the properties of metastable hadronic compact stars obtained 
within different relativistic mean-field models for the equation of state (EoS) of hadronic matter. In particular, we 
compare the predictions of the Quark-Meson Coupling (QMC) model [22|, [23| with those of the non-linear Walecka 
model (NLWM) [24j parametrizations given by Glendenning-Moszkowski (GM) (25| . 

For the quark phase we have adopted a phenomenological EOS [26J which is based on the MIT bag model for 
hadrons. The parameters here are: the mass m s of the strange quark, the so-called pressure of the vacuum B (bag 
constant) and the QCD structure constant a s . For all the quark matter model used in the present work, we take 
m u = nid = 0, m s = 150 MeV and a s = 0. 

In the QMC model quark degrees of freedom are explicitly taken into account: baryons are described as a system of 
non-overlapping MIT bags which interact through the effective scalar and vector mean fields. The coupling constants 
are defined at the quark level. An attractive aspect of the model is that different phases of hadronic matter, from very 
low to very high baryon densities, can be described within the same underlying model, namely the MIT bag model: 
matter at low densities is a system of nucleons interacting through meson fields, with quarks and gluons confined 
within MIT bag; at very high density one expects that baryons and mesons dissolve and the entire system of quarks 
and gluons becomes confined within a single, big MIT bag. 

In the case of the Glendenning-Moszkowski EoS (25[, we have paid special attention to the role played by the 
hyperon-meson couplings. In fact, all previous works on metastable hadronic stars [III, Ell EE Ell EH] have uniquely 
considered the case of "low" values for these quantities (x a — 0.6 for the ratio between the hyperon-u meson to 
nucleon-cr meson coupling) . As it is well known, larger values of the hyperon-meson couplings (constrained by the 
empirical binding energy of the A particle in nuclear matter) make the EoS stiffer and increase the value of the 
Oppenheimer-Volkoff mass for the hadronic stellar sequence [25j. In addition, as we demonstrate in the present work, 
increasing the values of the hyperon-meson couplings shifts the bulk transition point for quark deconfinement to 
higher densities and increments the value of the critical mass M cr (see ref. [H HIE! and Section Iffll for the explicit 
definition of this quantity) for the hadronic stellar sequence. Thus our study is relevant in connections with the recent 
measurements of highly-massive "neutron stars" mentioned above. 

A brief review of the NLW and QMC models is given in Section [TTJ The quantum nucleation of a quark matter 
drop inside hadronic matter is briefly reviewed in Section IlIII Our main results are presented in Section llVl whereas 
the main conclusions are given in Section [V] 

II. THE FORMALISM 

In the present section we review the models used in this work, namely the GM parametrizations [25| of the NLWM 
and the quark-meson coupling (QMC) model including hyperons. 

A. The non-linear Walecka model 

The Lagrangian density, including the baryonic octet, in terms of the scalar er, the vector-isoscalar w M and the 
vector-isovector meson fields reads (see e.g. 0,113, HI]) 
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£fiadrons H~~ ^leptons (1) 

where the hadronic contribution is 

^hadrons — ^baryons H~ ^mesons (2) 

with 

C baryons = $[-fD*-M%]1>, (3) 

baryons 

where 

D,j, = idp - guBU/j. - g P Bt~B ■ Pn, (4) 
and M B = Mb — g a B<3. The quantity ts designates the isospin of baryon B. The mesonic contribution reads 

^mesons £{X ~t~ ^uj ~\~ &pi (5) 

with 

= \(d^a - mla 2 ) + ^ + ±\a\ (6) 

= -^n fJl M ,1L ' + ^m 2 J uj t _ t uj fl , n MW = d^u v - d v u^ (7) 

£ P = -\b»v ■ + ]^n 2 p p^ -p^, B fiV = d fJ ,p v -d v p fl -g p (p fl x p„) (8) 
For the lepton contribution we take 

Cleptons = ^2 4>l (*7jU^ ~ m i'h (9) 

leptons 

where the sum is over electrons and muons. In uniform matter, we get for the baryon Fermi energy cfb — g^B^o + 



g P Bt3BP03 + y kp B + Mg , with the baryon effective mass M B = M — g a B<J- 

We will use the GM1 and GM3 parametrizations of NLWM [25| both fitted to the bulk properties of nuclear matter: 
for GM1 (GM3) the compressibility is 300 (240) MeV and the effective mass at saturation M* = 0.7 M (M* = 
0.78 M). The inclusion of hyperons involves new couplings, the hyperon-nucleon couplings: g a B = %oB 9a-, 9uB — 
Xu>b gu, 9pB = XpB 9p- For nucleons we take x a B, XuB, x p b — 1 and for hyperons we will consider the couplings 
proposed by Glendenning and Moszkowski [25j ] . They have considered the binding energy of the A in nuclear matter, 
B\, 

—j^j = ~ 28 MeV = ^ 9uj - x (j g a (J (10) 

to establish a relation between x a and x u . Moreover, known neutron star masses restrict x a to the range 0.6 — 0.8. 
We will take x p — x a and will consider x a = 0.6, 0.7, 0.8. 
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B. The quark- meson coupling model for hadronic matter 

In the QMC model, the nucleoli in nuclear medium is assumed to be a static spherical MIT bag in which quarks 
interact with the scalar and vector fields, a, to and p and these fields are treated as classical fields in the mean field 
approximation [12,[23|]. The quark field, ip q (x), inside the bag then satisfies the equation of motion: 



i p - (m° q - g%a) ~ g* u 7 ° + \g q p r z p w 



ip q (x) = , q = u,d,s, (11) 



where m° is the current quark mass, and g q , g^ and g q p and denote the quark- meson coupling constants. The 
normalized ground state for a quark in the bag is given by 

%{v, t) = N q cxp (-ie q t/R B ) ( . „ 1° [ Xqr / RB ) U A^=, (12) 
qy ' q y q/ ' \ ip q a ■ rji (x q r/R B ) J v ; 

where 

/ 1 \ lfl Q - R B m* 

e q = n q + R B [gl u, + -g q r zP03 j ; Pq = ^ +Rg ^ , (13) 

with the normalization factor given by 

A/ - " 2 = 2R%j 2 (x q ) [n q (n q - 1) + R B m* q /2] I x\ , (14) 



where VL q = y x q + (Rb ™j) 2 , rn* — mj — g q a, Rb is the bag radius of the baryon, and \ q is the quark spinor. The 

quantities ip q , e q , f3 q , J\f q , fl q , m* all depend on the baryon considered. The bag eigenvalue, x q , is determined by the 
boundary condition at the bag surface 

jo{x q ) = (3 g ji(x q ) . (15) 
The energy of a static bag describing baryon B consisting of three ground state quarks can be expressed as 

E h B ^ = J2n q ^-^ + -^R%B B> (16) 

where Z B is a parameter which accounts for zero-point motion and B B is the bag constant. The effective mass of a 
nucleon bag at rest is taken to be M B = E B &g . Th 
effective mass, M B with respect to the bag radius 



nucleon bag at rest is taken to be M B = E B &g . The equilibrium condition for the bag is obtained by minimizing the 



— # = . 17 

dR* B K ' 

For the QMC model, the equations of motion for the meson fields in uniform static matter are given by 

o = J2^BC B (a) 2 -^±l M b(°) k 2 dk (18) 

^ 2^ 2 J [k* + M B 2 (a)] 1/2 



ml 



B 



m 2 a, = J2 9^b (2J B + 1) 4/(6^ 2 ) , (19) 

B 



m 2 



Pos =Y,9 P BhB {2Jb + l)k B /(6n 2 ) . (20) 



In the above equations J Bl I 3B and k B are respectively the spin, isospin projection and the Fermi momentum of the 
baryon species B. For the hyperon couplings we take x w = 0.78 and x p — 0.7. The coupling x a is an output of the 
model and is approximately equal to 0.7. Note that the s-quark is unaffected by the a and u> mesons i.e. g# = = . 
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In Eq. JT8J) we have 

-,b. 



9aBC B {<y)= ^ = = 2^ n <i9lS B {<j) (21) 



where 



Aa S q U g {ftg - 1) + R B m*/2 

The total energy density and the pressure including the leptons can be obtained from the grand canonical potential 
and they read 



1 2 2 f 2 2 f 2 2 
£ = 2 m ° a + 2 m " w ° + 2 m P Po3 



£ 2J j^ I B k 2 dk [fc 2 + Mf(a)] 1/2 + £ 1 jf ' fc 2 dfc [fc 2 + m 2 ] 



21 V2 



(23) 



P 



1 2 2 , f 22,1 22 
: m„cr + -m,.,w + -m pPo3 



2 ° 2 

1 2 J B + 1 /"^ fc 4 dfc 1^ 1 /c 4 dfc 



3 ^ 2tt 2 Jo [ fc2 + m* 2 (ct)] 1/2 3 ; 71-2 Jo [jfca + m 2 ] 1/2 



(24) 



For the bag radius we take i?jy = 0.6 fm. The two unknowns Zn and P/v for nucleons are obtained by fitting the 
nucleon mass M = 939 MeV and enforcing the stability condition for the bag at free space. The values obtained are 
Z N = 3.98699 and B^ 4 = 211.303 MeV for m u = m d = MeV and Z N = 4.00506 and B X J 4 = 210.854 MeV for 
m u — md — 5.5 MeV. We take these bag values, B B , for all baryons and the parameter Z B and R B of the other 
baryons are obtained by reproducing their physical masses in free space and again enforcing the stability condition 
for their bags. Note that for a fixed bag value, the equilibrium condition in free space results in an increase of the 
bag radius and a decrease of the parameters Z B for the heavier baryons. The set of parameters used in the present 
work is given in Ref. [29| . 

Next we fit the quark-meson coupling constants g%, g^ = 3g^, and g p = g q for the nucleon to obtain the correct 
saturation properties of the nuclear matter, E B = e/p — M — —15.7 MeV at p — po — 0.15 fm~ 3 , a sym — 32.5 MeV, 
K = 257 MeV and M* = 0.774A/. We have g$ = 5.957, g uN = 8.981 and g pN = 8.651. We take the standard values 
for the meson masses, m a = 550 MeV, m w = 783 MeV m p = 770 MeV. 



III. QUANTUM NUCLEATION OF QUARK MATTER IN HADRON STARS 

Let us consider a pure hadronic star whose central pressure (density) is increasing due to spin-down or due to mass 
accretion (from a companion or from the interstellar medium) . As the central pressure approaches the deconfinement 
threshold pressure Po (see Fig. [2]), a drop of non /3-stable quark matter (Q*), but with flavor content equal to that 
of the /3-stable hadronic phase, can be formed in the central region of the star. The process of drop formation is 
regulated by its quantum fluctuations in the potential well created from the difference in the Gibbs free energies of 
the hadron and quark phases d, [ljl [H| 

U(TZ) = -7rn 6)Q . (n Q * - p H )TZ 3 + AircrK 2 (25) 

where 1Z is the radius of the Q* droplet (supposed to be spherical), n&,Q* is the quark baryon number density, pq* 
and p B are the quark and hadron chemical potentials at a fixed pressure P and a is the surface tension for the 
surface separating the hadron from the Q* phase. Notice that p is the same as the bulk Gibbs energy per baryon 
g = (P + e)/n B = (Y^ pin^)jn B . Notice also that we have neglected the term associated with the curvature energy, 
and also the terms connected with the electrostatic energy, since they are known to introduce small corrections 
[lol . \v\ . The value of the surface tension a for the interface separating the quark and hadron phase is poorly known, 
and typically values used in the literature range within 10 - 50 MeV far 2 00,^1. 
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The time needed to form the hrst drop (nucleation time) can be straightforwardly evaluated within a semi-classical 
approach [t| [nj. First one computes, in the Wentzel-Kramers-Brillouin (WKB) approximation, the ground state 
energy E and the oscillation frequency v§ of the drop in the potential well U(TZ). Then, the probability of tunneling 
is given by 



Po = exp 



A(Ep) 

h 



where A is the action under the potential barrier which in a relativistic framework reads 

A(E) = - [ J[2M(K)c 2 +E- U(K)][U(K) - E] , 
cJk- 



(26) 



(27) 



being 1Z± the classical turning points and 

M(K) = ^p H (l-^) 2 ll 3 (28) 
3 V n H J 

the droplet effective mass, with pn and uh the hadron energy density and the hadron baryon number density, 
respectively. The nucleation time is then equal to 

r - (voPoNc)- 1 , (29) 

where N c is the number of virtual centers of droplet formation in the star. A simple estimation gives iV c ~ 10 
0>[I3|- The uncertainty in the value of N c is expected to be within one or two orders of magnitude. In any case, all 
the qualitative features of our scenario will not be affected by this uncertainty. As a consequence of the surface effects 
it is necessary to have an overpressure AP = P — P > with respect to the bulk transition point Po to create a drop 
of deconfinement quark matter in the hadronic environment. The higher the overpressure, the easier to nucleate the 
first drop of Q* matter. In other words, the higher the mass of the metastable pure hadronic star, the shorter the 
time to nucleate a quark matter drop at the center of the star. 

In order to explore the astrophysical implications of quark matter nucleation, following ref. [Til [T3|. we introduce 
the concept of critical mass for the hadronic star sequence. The critical mass M cr is the value of the gravitational mass 
of a metastable hadronic star for which the nucleation time is equal to one year: M cr — Mhs{ t = lyr). Therefore, 
pure hadronic stars with Mhs > M cr are very unlikely to be observed, while pure hadronic stars with Mhs < M cr 
are safe with respect to a sudden transition to quark matter. Then M cr plays the role of an effective maximum mass 
for the hadronic branch of compact stars (see discussion in Ref. [12]). While the Oppenheimer-Volkov maximum mass 
is determined by the overall stiffness of the equation of state for hadronic matter, the value of M cr will depend in 
addition on the properties of the intermediate non /3-stable Q* phase. 



IV. RESULTS AND DISCUSSION 



In this section we present and discuss our results for stellar configurations obtained using the equation of state (EoS) 
models described in section II. In particular, we determine the region of the pure hadronic star sequence where these 
compact stars are metastable, the value of the corresponding critical mass M cr , and the final fate of this configuration 
after quark matter nucleation, i.e. whether it will evolve to a quark star or to a black hole. 

In Fig Q] the EoS for the models discussed are plotted for the range of densities of relevance for the discussion that 
follows. For GM1 and GM3 we have considered three different hyperon-meson coupling as discussed above. The 
QMC EoS corresponds approximately to x a — 0.7. A higher value of the hyperon couplings xi corresponds to stiffer 
EoSs: at high densities we have vector dominance defined by the magnitude of x u , x p . It is clear from Fig. [T]that 
the onset of hyperons (represented by the change of slope in the EoS curves) occurs for the smaller x a values at lower 
densities. The nucleonic EoS for QMC is very soft and therefore the onset of hyperons occurs at quite high densities, 
e = 373.87 MeV/fm 3 . As a consequence although QMC is softer than GM1 EoS at lower densities, it becomes, at 
higher densities, stiffer than GMl(a; CT = 0.6) and very close to GMl(i ff = 0.7). 

In Figs. [21 we plot the Gibbs' free energy per baryon for the hadronic phase and for the corresponding Q* phase 
using the various EoS models (couple of continuous and dashed curves with the same color) considered in the present 
work. It is clearly seen that in the case of the GM1 or GM3 EoS models, the lower the value of the hyperon coupling 
x a , the softer the EoS (see also Fig [l} and the lower the pressure Po at the crossing between the hadronic and the 
Q* phase. This will give rise to lower critical masses for the smaller x a values (see Tabl. I-III below). The Q* phase 
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FIG. 1: Hadronic Eos for QMC and for GM1 and GM3 with the different hyperon-meson couplings discussed in the text. 



is very sensitive to the particle content and it is due to this fact that, although in Fig. [2] the EoS for QMC is softer 
than the EoS for GM1 with x a — 0.8 and 0.7, its crossing with the Q* phase occurs at higher pressures. A similar 
observation occurs in the figure with the GM3 results. This behavior will reflect itself on values of the critical masses 
M cr . 

In Fig. [3J we show the mass-radius (MR) curve for pure HS within the QMC model for the EoS of the hadronic 
phase, and that for hybrid stars or strange stars for different values of the bag constant B. The configuration marked 
with an asterisk on the hadronic MR curves represents the hadronic star for which the central pressure is equal to 
the threshold value P an d the quark matter nucleation time is r = oo. The full circle on the hadronic star sequence 
represents the critical mass configuration, in the case a = 30 MeV/fm 2 . The full circle on the HyS (SS) mass-radius 
curve represents the hybrid (strange) star which is formed from the conversion of the hadronic star with Mhs = M cr . 
We assume |18| that during the stellar conversion process the total number of baryons in the star (or in other words 
the stellar baryonic mass) is conserved. Thus the total energy liberated in the stellar conversion is given by the 
difference between the gravitational mass of the initial hadronic star (Mj n = M cr ) and that of the final hybrid or 
strange stellar configuration with the same baryonic mass (Mf in = Mqs(M^ t ) ): 

E conv = (M m - M fln )c 2 . (30) 

As we can see from Fig. [3J for the case of the QMC model, the region of metastability of pure hadronic stars (the 
part of the MR curve between the asterisk and the full circle) is very narrow. For this hadronic EoS, the quark star 
sequence can be populated only in the case of "small" values of the bag constant (B < 80 MeV/fm 3 , in this case the 
final star is a strange star). In all the other cases the critical mass hadronic star will form a black hole. 

For comparison we plot the MR curve obtained with the GMI parametrization for the same surface tension (a — 
30MeV/fm 2 ) and two values for bag constant (B = 75 and 100 MeV/fm 3 ). We consider the two extreme values of 
the hyperon couplings studied in this work. The dots and stars have the same meaning as in Fig. [31 We see that 
for the cases plotted the only configuration that does not end in a black hole has the smallest bag constant and 
hyperon coupling considered. In the present model, however, the configuration with the central pressure Pq and the 
M cr configuration are quite separated, contrary to what was observed with QMC, Fig. [3) 

The larger mass difference between the star with the central pressure Po and the one with the M cr occurs when 
these stars have small masses. A small change in the central energy density corresponds to a large change in the 
mass. If instead of plotting the MR graph we would have plotted the corresponding mass-central pressure (MP) 
graph a larger difference between these two configurations would be expected. This is seen in Figs. [5] and [5] where 
the mass-pressure curves for the family of stars obtained respectively within QMC and GMI are plotted for two bag 
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FIG. 2: The Gibbs' energy per particle for the /3-stable hadronic phase (continuous curves) and for the respective Q* phase 
(dashed curves). The upper panel refers to the GM1 and the lower panel to the GM3 EoS. The results for the QMC model are 
plotted in both panels. 



constants and two values of the surface tension (a = 10 and 30 MeV/fm 2 ). We conclude that when the M cr star is 
almost on top of the Po star in the MR curves, these stars lie on or close to the plateau that contains the maximum 
mass configuration. A large separation between these two configurations corresponds to a phase transition which 
occurs during the rise of the MR curve before the plateau. Due to the softness of the QMC EOS, hyperons set on at 
quite large energy densities and the star with the central Po pressure only occurs at high densities. We also conclude 
that a smaller surface tension hastens the transition and the critical mass is closer to the Po mass. 

In Tables HI UTl and lllll we give the gravitational (M cr ) and baryonic (M°.) critical mass values for the hadronic star 
sequence, together with the central hyperon fraction (fy,cr = ny/ns, *.e. the ratio between the total hyperon number 
density and the total baryon number density at the center of the critical mass star). We also report the value of the 
gravitational mass (Mf in ) of the final quark star configuration and the total energy E conv = (M cr — Mf in )c 2 
released in the stellar conversion process, assuming baryon mass conservation (i.e. no matter ejection) [l8|. The 
gravitational (MQs, m ax) and the baryonic (Mq S max ) mass of the maximum mass configuration for the quark (hybrid 
or strange) star sequence are also included. The value of the latter quantity is relevant to establish whether the 
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FIG. 3: Mass-radius relation for a pure HS described within the QMC model and that of the HyS or SS configurations for 
several values of the bag constant and m 3 = 150 MeV and a a — 0. The configuration marked with an asterisk represents in all 
cases the HS for which the central pressure is equal to Po- The conversion process of the HS, with a gravitational mass equal 
to M cr , into a final HyS or SS is denoted by the full circles connected by an arrow. In all the panels a is taken equal to 30 
MeV/fm 2 . 



critical mass hadronic star will evolve to a quark star {M h cr < MQ S max ) or will form a black hole (M h cr > MQ S max ). 
The entries in Tables HI lU [49] . and lllll are relative respectively to the GM1, GM3 and QMC equation of state for the 
hadronic phase. For the quark phase we consider four different values of bag constants, 75, 85, 100 and 150 MeV/fm 3 , 
and two different values for quark-hadron surface tension, 10 and 30 MeV. Notice that for the quark matter parameter 
set adopted in the present work (see Section I), strange quark matter is absolutely stable [H, [34[ only for B = 75 
MeV/fm 3 . 

Some comments are in order: the critical masses increase with the increase of the hyperon couplings. This increase 
can be as large as 0.3 - 0.4 M Q when x„ changes from 0.6 to 0.8; the critical mass is also dependent on the particle 
content, namely of the strangeness content, and this explains the different relative positions for the different bag 
pressures of the QMC result which essentially corresponds to x„ = 0.7. Due to the fact that the EoS for QMC is 
very soft, the hyperon onset occurs at quite high densities and therefore the critical mass is always quite high for 
this model. The critical mass increases with the bag constant because a larger bag constant corresponds to a stiffer 
quark EoS and therefore the phase transition to the quark phase will occur at larger densities. When the critical mass 
hadronic star is converted to a black hole, this is indicated in Tables HI HT1 and IIIII with a entry BH, in the columns 
for Mfi n and E conv (no energy will be radiated as soon as the star pass the event horizon). Notice that, in the case 
of the GM3 model with x a — 0.6 and B = 150 MeV/fm 3 there is no entry for the critical mass value (and for M% r , 
Mfi n and E conv ) since in this case the nucleation time of the maximum mass hadronic star (MHS,max) is much larger 
than one year (i.e. the star is metastable with a life-time comparable or much higher than the age of the universe). 

We observe from the results in tables Q] UH that increasing the value of the hyperon coupling constants (for fixed B 
and a) reduces the central hyperon fraction (fy,cr) of the critical mass star, and increases the energy released during 
the conversion into a quark or hybrid star (for those configurations which will not form a black hole). 

In Fig. we show the internal composition for the hadronic star with a gravitational mass M = 2.081M Q and 
radius R = 12.6 km, obtained using the GM1 parametrization with x a — 0.8. This star corresponds to the critical 
mass configuration when we consider B = 150 MeV/fm 3 and a = 30 MeV/fm 2 (see table HJ). As we see, this star has a 
considerable central hyperon fraction (/y, cr = 0.299) and a wide hyperonic matter core which extend up to Ry ~ 8.7 
km. On the top of this core, one has a nuclear matter layer (Ry < r < R cr ust) with a thickness of about 3.4 km. The 
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FIG. 4: Mass-radius relation for a pure HS described within the GM1 parametrization and that of the HyS or SS configurations 
for two values of the bag constant (B — 75 and 100 MeV/fm 3 ) and two values of the hyperon-meson coupling [x a — 0.6, and 
0.8) and m s = 150 MeV and a s — 0. The configuration marked with an asterisk represents in all cases the HS for which the 
central pressure is equal to Po- The conversion process of the HS, with a gravitational mass equal to M cr , into a final HyS or 
SS is denoted by the full circles connected by an arrow. In all the panels a is taken equal to 30 MeV/fm 2 . 



stellar crust extends from R crU st up to R. 

It has been argued by several authors [HI, H(| H3, [38| that if strange quark matter (SQM) is absolutely stable 
[33|,[34|], then all compact stars are likely to be strange stars. The argument in favor of this thesis is the following: 
if the interstellar medium is sufficiently contaminated by quark nuggets (i.e. lumps of SQM), then the presence of a 
single quark nugget in the interior of a "normal" neutron star (hadronic star) is sufficient to trigger the conversion 
of the star to a strange star [34], |35| . Likely the quark nug get contamination of the interstellar medium is the result 
of the merging of strange stars in binary systems [3a, l39j |. Under these conditions, compact star progenitors could 
capture a quark nugget during their lives (i.e. during the various nuclear burning stages of the stellar evolution). 
Thus, according to this argument, a strange quark seed will be present in all new born compact stars, and thus the 
conversion to a strange star will happen immediately, without a metastable hadronic star being formed first. This 
is a plausible scenario, however it would be relevant only for a few of the stellar models considered in our work, 
i.e. those relative to the value B = 75 MeV/fm 3 for the bag constant (see Tables HI |TT1 and IIII[) for which SQM is 
absolutely stable, and will not have any effect upon the existence of metastable hadronic compact stars in all the 
other cases considered in the present work. The magnitude of the flux of q uark nuggets in the interstellar medium 
(which is a crucial quantity for the validity of the scenario of ref.s HflHfilHJHl]) nas ^ een estimated (36| making the 
assumption that all pulsars exhibiting glitches must be "normal" neutron stars (hadronic stars), not strange stars. 
This assumption is based primarily on the nearly total lack of models for the glitch phenomenon with strange stars 
(see anyhow ref. EftS] ), while such models have been quite successfully developed in the case of hadronic stars (see 
e.g. ref. [H, El, Hj, [45 1 ) . However, recent studies have established the possibility of an inhomogeneous crystalline 
color superconducting phase (LOFF phase) in the interior of strange stars (see [46],[47[ and references therein quoted), 
or the likely existence of a SQM crystalline crust in strange stars [48[. These late theoretical developments rise the 
possibility to explain pulsar glitches with strange star based models, and thus require as useful a recalculation of the 
astrophysical limits of the flux of quark nuggets. The scenario discussed in the present work is an alternative to the 
scenario [35l.l36l.l37j. 38] according to which all compact stars are strange stars, which requires that SQM is absolutely 
stable. 
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FIG. 5: Mass-pressure relation for a pure HS described within the QMC model and that of the HyS or SS configurations 
two values of the Bag constant(75 and 100 MeV/fm 3 ) and m s = 150 MeV and a s — 0. The configuration marked with an 
asterisk represents in all cases the HS for which the central pressure is equal to Po- The conversion process of the HS, with a 
gravitational mass equal to M cr , into a final HyS or SS is denoted by the full circles connected by an arrow. Two values of the 
surface energy a were considered 10 MeV/fm 2 (left) and 30 MeV/fm 2 (right). 



V. CONCLUSIONS 

It has been recently shown [ll|, [l^, [IH HH, EH that pure hadronic compact stars, above a threshold value of their 
gravitational mass, are metastable to the conversion to quark stars. In this work we have done a systematic study of 
the metastability of pure hadronic compact stars using different relativistic hadronic models for the equation of state 
of hadronic dense matter. In particular, we have used and compared the quark-meson coupling (QMC) model with 
those for the Glendenning-Moszkowski parametrization of the non-linear Walecka model (NLWM). In the case of the 
QMC model, we have obtained that the region of metastability of pure hadronic stars is very narrow. For the GM 
model, we have investigated the effect of the hyperon couplings on the critical mass of the hadronic star sequence and 
on the stellar conversion energy. We have found that increasing the value of the hyperon coupling constants shifts 
the bulk transition point for quark deconfincmcnt to higher densities, increasing the value of the critical mass for the 
hadronic stellar sequence, and thus makes the formation of quark stars less likely. The nucleonic EoS for QMC is very 
soft and therefore the onset of hyperons occurs at quite high densities, which gives rise to large critical masses. The 
conversion to a quark star will occur only for a small value of the bag constant. Finally we point out that both QMC 
and GM1 with the largest values of the hyperon-meson couplings predict limiting masses [12j | which may be as high 
as 1.9 - 2.1 M Q . 

These values would be able to describe highly-massive compact stars, such as the one associated to the millisecond 
pulsars PSR B1516+02B HI, and nearly the one in PSR J1748-2021B [2ljj . 
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FIG. 6: The same as Fig [5] for the GM1 parametrization with the x a — 0.6 hyperon coupling. 
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TABLE I: Critical masses and energy released in the conversion process of an HS into a QS, for several values of the hyperon 
coupling x a , of the Bag constant B and the surface tension a. The GM1 parameter set has been used for the hadronic EoS. 
Column labeled MQs,max {MQ Smax ) denotes the maximum gravitational (baryonic) mass of the final QS sequence. The value 
of the critical gravitational (baryonic) mass of the initial HS is reported on column labeled M cr (M^ r ) whereas those of the 
mass of the final QS and the energy released in the stellar conversion process are shown on columns labeled M/,„ and E conv 
respectively. BH denotes those cases in which due to the conversion the initial HS collapses into a black hole. fy,cr = ny/ns 
denotes central hyperon fraction of the critical mass star (i.e. the ratio between the total hyperon number density and the total 
baryon number density at the center of the critical mass star). Units of B and a are MeV/fm 3 and MeV/fm 2 respectively. All 
masses are given in solar mass units and the energy released is given in units of 10 51 erg. m s and a s are always taken equal to 
150 MeV and respectively. The Oppenheimer-Volkoff maximum masses for pure hadronic stars in the case of the GM1 EoS 
are: M H s,max = 1.790M S (when x a = 0.6), 1.996 M (x a = 0.7), and 2.169 M Q {x a = 0.8). 



a = 10 a = 30 





B 


MQS,max 




M cr 


Ml f Y ,cr 


M fin 


R 


M cr 


Ml fy,cr 


M fin 


J^conv 


0.6 


75 


1.630 


1.968 


1.326 


1.454 0.079 


1.254 


128.4 


1.471 


1.630 0.147 


1.387 


149.4 




85 


1.542 


1.812 


1.447 1.596 0.134 1.385 


110.4 


1.540 


1.711 0.201 


1.479 


125.5 




100 


1.457 


1.661 


1.598 


1.789 0.261 


BH 


BH 


1.658 


1.865 0.332 


BH 


BH 




150 


1.447 


1.601 


1.770 


2.010 0.527 


BH 


BH 


1.790 


2.036 0.637 


BH 


BH 


0.7 


75 


1.630 


1.968 


1.442 


1.595 0.059 


1.361 


144.8 


1.602 


1.794 0.119 


1.507 


169.6 




85 


1.542 


1.812 


1.584 


1.764 0.111 


1.509 


134.2 


1.685 


1.892 0.170 


BH 


BH 




100 


1.457 


1.661 


1.724 


1.950 0.197 


BH 


BH 


1.791 


2.036 0.253 


BH 


BH 




150 


1.518 


1.686 


1.905 2.188 0.370 


BH 


BH 


1.931 


2.223 0.403 


BH 


BH 


0.8 


75 


1.630 


1.968 


1.592 


1.782 0.044 


1.498 


167.9 


1.763 


2.000 0.095 


BH 


BH 




85 


1.542 


1.812 


1.735 


1.954 0.085 


BH 


BH 


1.841 


2.091 0.132 


BH 


BH 




100 


1.457 


1.661 


1.879 


2.152 0.152 


BH 


BH 


1.946 


2.243 0.194 


BH 


BH 




150 


1.518 


1.686 


2.054 


2.391 0.275 


BH 


BH 


2.081 


2.429 0.299 


BH 


BH 



15 



TABLE II: Same as Tab. I, but for the GM3 parameter set for the hadronic EoS. The Oppenheimer-Volkoff maximum masses 
for pure hadronic stars in the case of the GM3 EoS are: M H s,max = 1.554M Q (when x a = 0.6), 1.732 Mq (x a = 0.7), and 
1.875 M (av = 0.8). 



a = 10 a = 30 





B 


M QS , max 






M b cr f Y , cr 


M fin 


E conv 


M cr 


Ml fy, cr 


M fin 


Econv 


0.6 


75 


1.630 


1.968 


1.237 


1.351 0.092 


1.175 


111.6 


1.269 


1.389 0.110 


1.204 


115.4 




85 


1.543 


1.812 


1.350 


1.482 0.165 


1.298 


91.3 


1.362 


1.497 0.178 


1.310 


92.8 




100 


1.465 


1.673 


1.461 


1.626 0.307 


1.431 


54.9 


1.469 


1.636 0.316 


1.438 


55.8 




150 


1.487 


1.658 


















0.7 


75 


1.630 


1.968 


1.373 


1.510 0.078 


1.297 


136.4 


1.402 


1.545 0.091 


1.324 


140.4 




85 


1.543 


1.812 


1.511 


1.680 0.157 


1.447 


113.2 


1.541 


1.717 0.182 


1.475 


117.6 




100 


1.465 


1.673 


1.610 


1.806 0.253 


BH 


BH 


1.645 


1.851 0.295 


BH 


BH 




150 


1.495 


1.667 


1.716 


1.945 0.419 


BH 


BH 


1.723 


1.956 0.444 


BH 


BH 


0.8 


75 


1.630 


1.968 


1.574 


1.759 0.076 


1.482 


165.2 


1.611 


1.806 0.092 


1.516 


170.7 




85 


1.543 


1.812 


1.694 


1.913 0.143 


BH 


BH 


1.744 


1.979 0.181 


BH 


BH 




100 


1.465 


1.673 


1.771 


2.014 0.204 


BH 


BH 


1.802 


2.057 0.234 


BH 


BH 




150 


1.495 


1.668 


1.848 2.119 0.295 


BH 


BH 


1.856 


2.131 0.309 


BH 


BH 



TABLE III: Same as Tab. I, but for the QMC parameter set for the hadronic EoS. The Oppenheimer-Volkoff maximum mass 
for pure hadronic stars in the case of the QMC EoS is: MHS,max = 1.927M©. 



a = 10 a = 30 



B 


MQS,max 


MQS,max 


M cr Ml f Y ,cr 


M fin 




M cr 


Ml 


fY,cr 


M fm 


TP 


75 


1.630 


1.968 


1.587 1.768 0.044 


1.488 


176.8 


1.694 


1.903 


0.090 


1.585 


195.8 


85 


1.530 


1.793 


1.705 1.917 0.096 


BH 


BH 


1.768 


1.998 


0.145 


BH 


BH 


100 


1.454 


1.656 


1.790 2.027 0.168 


BH 


BH 


1.830 


2.080 


0.211 


BH 


BH 


150 


1.479 


1.638 


1.898 2.171 0.352 


BH 


BH 


1.909 


2.187 0.377 


BH 


BH 



